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Abstract 

Neutrino telescopes of large area offer the possibility of searching for indirect 
signals of relic neutralinos in the galactic halo, due to annihilations in the 
Sun or the Earth. Here we investigate the sensitivity, using a supergravity 
scheme where the soft scalar mass terms are not constrained by universality 
conditions at the grand unification scale. We first discuss in which regions of 
the supersymmetric parameter space the neutralino may be considered as a 
good candidate for cold dark matter. The discovery potential of the search 
using neutrino telescopes is then compared to that of the direct search for 
relic neutralinos. 



I. INTRODUCTION 

In large regions of the supersymmetric parameter space the neutralino turns out to be 
the Lightest Supersymmetric Particle (LSP) and, as such, it is stable, provided R-parity is 
conserved. Under these hypotheses the neutralino would have decoupled from the initial 
plasma in the early stages of the Universe, and would now be present as a relic particle [|TJ . 
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The LSP is a candidate for the Cold Dark Matter (CDM) that is believed to have played 
a key role in the formation of structures in the Universe, such as galaxies and clusters. 
Although there are other candidates for CDM, such as the axion or the axino, and scenarios 
for structure formation that do not involve a large density of CDM, the possibility that the 
neutralino provides most of the CDM remains in our opinion the most attractive option. 
The "Standard Model" of structure formation used to be one with an initially flat Harrison- 
Zeldovich spectrum of inflationary fluctuations and (essentially) the critical density of CDM: 
VLcdm = Pcdm/ Pcrit — 1- However, the advent of COBE and other data have suggested 
that this model needs to be modified, and there are three main contenders on the market 

One is a "Mixed" Dark Matter (or Cold Hot Dark Matter) model (CHDM), in which 
VLcdm — 0.7, there is a Hot Dark Matter component with VLhdm — 0.2, and baryons 
contribute VIb < 0.1, as reviewed in Section 2. Another is a model (ACDM) with a significant 
cosmological constant ~ 0.7 and Qcdm — 0-3. Finally, we mention a Cold Dark Matter 
with a tilted spectrum of initial fluctuations (TCDM), in which Qcdm — 1 is still possible, 
and an open model with Qcdm — 0.3. Calculations of the relic LSP abundance actually 
lead to values for the product Q^sph 2 , where h is the present Hubble expansion rate Hq in 
units of 100 Km - s _1 ■ Mpc -1 . There is still some observational uncertainty in this quantity, 
which may lie in the range 0.5 < h < 0.9. Within this range, the age of the Universe favours 
smaller values of h in the CHDM and TCDM scenarios, whereas larger values are possible in 
the ACDM scenario. Combining the estimates of Qcdm and h in each of the three scenarios, 
we find the preferred range 

QcDNih 2 = 0.2 ±0.1 (1) 

However, we repeat that not all the CDM need be constituted of LSPs, so Qlsp could in 
principle lie below the range (JJ). 

Even if relic neutralinos do not provide a significant fraction of dark matter, experimental 
evidence for them would add a relevant new piece of information on the early stages of the 
Universe. Various strategies for detection of relic neutralinos are currently being pursued || . 
The most straightforward technique (direct detection) consists in measuring the effect that 
an impinging neutralino may produce in an appropriate detector by its elastic scattering off 
a target nucleus [|J. Among the indirect ways of detecting relic neutralinos, one of the most 
promising ones is the observation, using neutrino telescopes of large area, of the up-going 
muons which would be generated by neutrinos produced by pair annihilation of neutralinos 
captured and accumulated inside celestial bodies such as the Earth and the Sun |^[7|]. 

At present, theory is unfortunately unable to offer firm predictions for the event rates for 
detection of relic neutralinos, since supersymmetric theories are still awaiting experimental 
verification. Only some hints for possible supersymmetric effects are available from acceler- 
ator data: supersymmetric theories would favor the unification feature of the gauge running 
constants || and the apparently relatively light Higgs boson mass @. Apart from these 
properties, the physical (correct) supersymmetric scheme is not known yet, and thus one 
has to consider a number of various possible scenarios. The detection rates for neutralinos 
depend very sensitively on the different supersymmetric schemes employed in the analysis. 

The least constrained theoretical model is represented by the Minimal Supersymmetric 
extension of the Standard Model (MSSM), which incorporates the same gauge group as the 
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Standard Model and the minimal supersymmetric extension of its particle content ||10|| . The 
Higgs sector contains two doublets H%, H 2 which give masses to down- and up- type quarks, 
respectively. This scheme provides a very useful framework for analyzing the phenomenol- 
ogy at the Mz scale with a minimal number of model-dependent restrictions. The main 
inconvenience of this approach is that one typically has to deal with a large number of free 
parameters. 

Much more ambitious are theoretical schemes where features at the Mz scale are derived 
from properties at the Grand Unification (GU) scale (Mqut), the link being provided by 
the Renormalization Group Equations (RGE's). One of the most attractive supersymmetric 
models is the one in which Electro- Weak Symmetry Breaking (EWSB) is induced radiatively 
| Till . This model is in no way the only possible model, but has the very nice feature of 



connecting the EWSB to soft supersymmetry breaking. Furthermore, the requirement of 
radiative EWSB is effective in reducing the number of the free parameters. 

In order to constrain further the scheme, one often makes a number of rather restrictive 
hypotheses. Typically one assumes that not only the gauge couplings but also the Yukawa 
couplings of b and r and the soft-breaking mass parameters (gaugino masses, scalar masses 
and trilinear couplings) unify at a GUT scale Mqut = O(10 16 GeV). These assumptions 
entail very strong consequences for neutralino phenomenology, and in particular for the 
properties of neutralino dark matter. However, it has been shown that many aspects of 
neutralino phenomenology may change quite significantly, if one relaxes the universality 



requirements [|T^,|T^]. This possibility has been explored in Ref. [J4| in an analysis of the 
neutralino relic abundance and of the event rates for direct detection in a wide range of the 
supersymmetric parameter space. 

In the present paper we extend the analysis of Ref. [14] to the evaluation of the event rates 
for indirect relic neutralino searches using neutrino telescopes. Our results are compared 
with new, more stringent experimental bounds, obtained using the Baksan detector |15| . 



Furthermore, we also present a comparison between the discovery potential of searches using 
neutrino telescopes and the direct method. 

The plan of the paper is as follows. In Sect. 2 we expand the above discussion of the 
CDM density, and derive the estimate ([!]) for fl C DMh 2 , by considering various cosmological 
models. The main features of the supersymmetric scheme employed in this paper are briefly 
described in Sect. 3, and in Sect. 4 we present our results on the neutralino relic abundance. 
Sect. 5 is devoted to the derivation of the flux of up-going muons due to neutralino- neutralino 
annihilation in the Earth and in the Sun. In Sect. 5 we also discuss some relevant features 
of the experimental layout required for the detection of the signals under study. Finally, 
results and conclusions are given in Sect. 6. 



II. DENSITY OF CDM IN COSMOLOGICAL MODELS 

Here we consider the neutralino as a CDM particle in the framework of several cosmo- 
logical models. Such models are primarily characterized by two dimensionless parameters 
h = F o /(100 Km s" 1 Mpc" 1 ) and Q = p/p c , where p is the relic cosmological density and 
p c ~ 1.88 ■ I0~ 29 h 2 g/cm 3 is the critical density. 
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Different measurements of the Hubble constant imply 0.4 < h < 1 [16|, but the re- 



cent measurements of extragalactic Cepheids in the Virgo |17| and Coma [18| clusters have 



narrowed this interval to 0.5 < h < 0.9. However, this range should be taken with some 
caution, because of the uncertainties involved in these difficult measurements. In particular, 
the value h = 0.5 has to be considered as only marginally allowed. Inspired mostly by theo- 
retical motivations (flatness problem and the beauty of the inflationary scenario), $7 = 1 is 
usually assumed. This value is consistent with the IRAS [TP]] data and the POTENT 
analysis, and no observational data contradict this value significantly. 

Dark Matter can be subdivided into baryonic DM, hot DM (HDM) and CDM. The den- 
sity of baryonic matter found from nucleosynthesis is taken [21J as Qh 2 = 0.025 ± 0.005. 



The other DM components are defined as hot or cold components depending on their veloc- 
ities at the moment when galaxies cross the horizon scale. If particles are relativistic they 
are called HDM particles, if not, CDM particles. The natural candidate for HDM is the 
heaviest neutrino, most probably the r neutrino. Structure formation in the Universe puts 
strong restrictions to the properties of DM in the Universe. A Universe with only HDM and 
baryonic DM gives a wrong prediction for the spectrum of fluctuations as compared with 



the measurements of COBE g2[, IRAS |L9| and the CfA jEJ survey. CDM and baryonic 
matter alone may explain the spectrum of fluctuations if the total density Q ~ 0.3. There 
is one more possible form of energy density in the Universe, namely the vacuum energy 
described by the cosmological constant A. The corresponding energy density is given by 
= A/(3Hq). Quasar lensing and the COBE results restrict the vacuum energy density 
to n A < 0.7. 

There are several cosmological models based on the four types of DM described above 
(baryonic DM, HDM, CDM and vacuum energy). These models predict different spectra of 
fluctuations to be compared with the data of COBE, IRAS, the CfA survey, etc. They also 
produce different cluster-cluster correlations, number densities of clusters, velocity disper- 
sions and other properties. The simplest and most attractive model for a correct description 
of all these phenomena is the mixed model (CHDM). This model is characterized by the 
following parameters: 

= 0, fi = fife + &CDAI + &HDM = 1) 

H ~ 50 Km s" 1 Mpc -1 (/i ~ 0.5), 
Vcdm ■■ Mhdm ■■ V b ~ 0.7 : 0.2 : 0.1, (2) 



where Qhdm — 0.2 is obtained in Ref. [24| from data on damped Lyman a clouds. In this 



CHDM model the central value for the CDM density is given by 

n CD Mh 2 = 0.18 (3) 

with uncertainties which may be estimated as < ±0.1. As already mentioned, the best 
candidate for the HDM particle is the r neutrino, and in the CHDM model with Q u = 0.2 
its mass should be m Vr ~ 4.7 eV. A very good fit to the cosmological data is given by a 
CHDM model where the HDM is constituted of two neutrinos {Cv 2 DM [^,^]). In our 
view, in either case, the most plausible candidate for the CDM particle is probably the 
neutralino (x). 
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In the light of recent measurements of the Hubble constant, the CHDM model faces 
a possible age problem. A lower limit on the age of the Universe to > 13 Gyr from 
globular clusters imposes an upper limit on the Hubble constant in the CHDM model: 
H < 50 Km s _1 Mpc -1 . This value is in slight contradiction with the recent observations 
of extragalactic Cepheids, which can be summarized as H > 60 Km s _1 Mpc -1 . However, 
it is too early to consider this as a serious conflict J27|| , if we take into account the many 
uncertainties and physical possibilities (e.g., the Universe can be locally overdense - see the 
discussion in ref. f26j). 

The age problem, if taken seriously, can be solved with the help of another successful 
cosmological model, ACDM. This model assumes that Q — 1 is provided by the vacuum 
energy (cosmological constant A) and CDM. In this case we have Qa — 0.7, Qcdm — 0.3 
and h < 0.7. Thus this model predicts VLcDuh 2 — 0.15 with an uncertainty of order 0.1. 

Finally, we mention two other CDM models. The first one is the tilted CDM model 
(TCDM), where the initial fluctuation spectrum is steeper than the Harrison- Zeldovich 
spectrum, and the second one is the CDM model with Q = Qcdm = 0.3 (CDM). Both 
models give a good fit to the observed spectrum of fluctuations as well as good agreement 
with the cluster data, though they may conflict with conventional inflationary prejudices. 

In Table 1 we summarize the estimates of flcDM for all these models and the ensuing 
values for QcDMh 2 . Taking into account the uncertainties, we conclude as mentioned in 

TABLE I. ^IcDhdh 2 for five cosmological models 





&-CDM 


h 


^CDhlh 2 


CHDM 


0.7 


0.5 


0.18 


ACDM 


0.3 


0.7-0.8 


0.15-0.19 


Cv 2 DM 


0.7 


0.5 


0.18 


TCDM 


1.0 


0.5 


0.25 


CDM 


0.3 


0.7-0.8 


0.15-0.19 



Sect. 1 that, for all the cosmological models considered here 

VcDMh 2 = 0.2 ±0.1 (4) 

In the following we will emphasize regions of the supersymmetric parameter space which 
yield a neutralino relic abundance Q x h 2 within the range of Eq.(|j), but considering also 
regions with lower Q x h 2 , since there can be additional forms of CDM, such as axions. 

III. THEORETICAL FRAMEWORK 

We turn now to a short presentation of the theoretical model employed here to describe 
the neutralino. We adopt a supersymmetric model whose essential elements are provided 
by a Yang-Mills Lagrangian, the superpotential, which contains all the Yukawa interactions 
between the standard and supersymmetric fields, and by a soft-breaking Lagrangian with 
the usual trilinear couplings (with parameters A's), the Higgs-mixing term /i, and the mass 
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terms (M, for the gaugino masses and rrii for the scalar masses). While, for simplicity, 
unification conditions at M GUT are imposed on gaugino masses: Mi(M GUT ) = m 1 / 2 , and 
for the trilinear couplings (A Q being their common value at Mqut), soft scalar masses are 
allowed to deviate from strict universality. More specifically, we consider here a departure 
from universality in the scalar masses at Mqut that splits the soft-supersymmetry-breaking 
mass parameters of the two Higgs doublets Mn x , M# 2 in the following way 

M 2 Hi (M GUT ) = m 2 (l + 5 t ) . (5) 

The parameters 5j are varied in the range (—1, +1), but are taken to be independent of the 
other supersymmetric parameters. 

Our supersymmetric parameter space is then constrained by a number of conditions: a) 
all experimental bounds on Higgs, neutralino, chargino and sfermion masses are satisfied 
(taking into account also the new data from LEP1.5 [p8j1 , b) the neutralino is the Lightest 
Supersymmetric Particle (LSP), c) constraints on the b — > s'j process and on the mass of 
the bottom quark m& assuming b-r Yukawa unification are satisfied, d) EWSB is realized 
radiatively, e) radiative EWSB occurs without excessive fine-tuning, f) the neutralino relic 
abundance does not exceed the cosmological bound. In particular, the requirements of 
radiative EWSB and of the universality conditions on the gaugino masses and on the trilinear 
couplings allow a reduction of the independent parameters to the following set (apart from 
the Si's): rai/2, rrio, A , tan/3 (tan/5 is the ratio v 2 /vi of the vacuum expectation values 
of the two Higgs doublets). It has to be emphasized that, because of the assumption of 
radiative EWSB, the parameter /i of the Higgs-mixing term in the superpotential is not a 
further independent parameter, but is a function of the previous set of parameters. 



As we emphasized in our previous paper [14], attention must be paid to the condition 



that radiative EWSB is satisfied without excessive tuning pU| , |5(]| . In Ref . [|TJJ] we found that 
requiring accidental cancellations among various competing terms not to exceed the 1% level 
sets the bound m x < 200 GeV. Further details of the theoretical scheme adopted here can 



be found in Ref. 14 



A departure from m universality of the type given in Eq. fl5|) may modify the neutralino 
phenomenology in a significant way. Two key parameters whose values may change sizably, 
as functions of the <Vs, are fi and the mass rriA of the CP-odd Higgs boson A. In turn, 
variations in fi and Ma may induce significant modifications in the neutralino properties. It 



is convenient to express these two parameters in the following way |L4] 



M 2 

/x 2 = J x m\ j2 + J 2 m\ + J 3 Alm 2 + J 4 A m m 1/2 — (6) 



and 



M\ = K x m\ l2 + K 2 ml + K 3 A 2 m 2 + K 4 A m m 1/2 - M 2 Z , (7) 



where the coefficients Ji and Ki, which are functions of tan/3 and of the 5i (except for J\ and 
K\ which depend on tan/3 only), are obtained from the RGE's. The coefficients J\, J 2 and 
Ki, K 2 are given in Fig. 1 for the case of m universality. From now on we will set A = 0, 
thus the other coefficients are irrelevant for our discussion. Except for very small values of 
tan/3 (tan/3 < 4), one has J 2 <^ Ji, which in turn implies a strong mi/ 2 -fi correlation and 
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a gaugino-like neutralino. However, even moderate departures from mo universality may 



modify this picture [p!2|,|14| . For instance, at large tan/3, non vanishing values of the 5,'s may 
yield a sizeable value of \J 2 \, which in turn entails either a more pronounced gaugino-like 
neutralino, when J 2 > 0, or a mixed higgsino-gaugino composition for the neutralino, when 
J 2 < 0. 

Also the coefficient K 2 plays a key role in establishing important phenomenological prop- 
erties. In the universal case (see Fig.l), K 2 is positive and sizeable, except at very large 
tan/? (i.e., tan/3 ~ 50). Thus, M4 turns out to be large except at very large values of tan/3, 
where it may approach the present experimental bound Ma > 55 GeV. Again, deviations 
from m universality may modify K 2 substantially and thus may change the value of Ma, 
too. We will see in the following how these properties affect the values of some important 
quantities, such as the relic abundance. 

For further details about the procedure we adopted to solve the RGE's and to implement 



the constraints due to b — > 57 and to m b , we refer to our previous paper |TJ]]. Here we 



only recall that we used 1-loop beta functions including the whole supersymmetric particle 
spectrum from the GUT scale down to Mz, neglecting the possible effects of intermediate 
thresholds. Two-loop and threshold effects, which are known to be of key importance 
for specific refined calculations, such as the unification of the gauge couplings, were not 
included, since we are interested in overall neutralino properties studied over a wide range 
of variation for the high-scale parameters. We have allowed generous ranges for b —>■ sj and 
rrib, to accommodate uncertainties in QCD corrections and the correct Yukawa unification 
condition, respectively. Furthermore, some physical solutions are unstable with respect to 
allowable variations in the strong gauge coupling constant a s . Although the best way of 
proceeding would be to allow a s to vary over its whole physical range, here, for simplicity, 
we have preferred to show our results only for some representative values of a s . In one case 
in Sect. 4 we show comparatively our results for two different values of the strong coupling 
constant. 

Our results for Q x h 2 and for the detection rates have been obtained by varying the 
parameters mo and vn,\i 2 on an equally-spaced linear grid over the ranges 10 GeV < < 
2 TeV, 45 GeV < m\i 2 < 500 GeV at fixed tan/3, with A Q = 0. Furthermore, we remark 
that all evaluations presented here are for positive values of /x, since negative values of \i are 
disfavoured by the constraints due to m& and the b — > 37 process. 

As far as the values of tan/3 are concerned, we note that recent global fits of the elec- 
troweak data within the MSSM |H| focused interest on two narrow intervals for tan/3: (i) 
very small values, tan/3 ~ 1 — 2 (i.e., close to the quasi-infrared fixed point for a given value 
of m t ), or (ii) very large values of tan/3 (i.e., of order mt/m^). In the present paper, as 
representative values for case (i) and case (ii), we take the values tan/? = 1.5 and tan/3 = 53, 
respectively. 



We conclude this section by summarizing the basic features of our model (i) the 
universality condition for the scalar masses is relaxed, (ii) rather relaxed restrictions from 
b — *■ 57 and the mass of the b quark are allowed, (iii) the fine-tuning condition limits the 



neutralino mass to m Y < 200 GeV. 



7 



IV. NEUTRALINO RELIC ABUNDANCE 



For the evaluation of Q x h 2 we have followed the standard method f|32]-|35||. We re- 



call that Vt x h 2 is essentially given by Vt x h 2 oc< cr ^ v > 1 : n J, where < a ann v >j nt is the 
thermally-averaged annihilation cross section, integrated from the freeze-out temperature 
to the present temperature. Then the key quantity to be evaluated is the annihilation 
cross-section. 

In the evaluation of cr ann we have considered the following set of final states: (1) fermion- 
antifermion pairs, (2) pairs of charged Higgs bosons, (3) one Higgs boson and one gauge 
boson, (5) pairs of gauge bosons. For the final state (1), the following diagrams have been 
considered: Higgs- and Z-exchange diagrams in the s channel and / exchange in the t 
channel. For the final states (2 to 5) we have included Higgs-exchange and Z-exchange 
diagrams in the s channel, and either neutralino (the full set of the four mass eigenstates) or 
chargino exchange in the t channel, depending on the electric charges of the final particles 

The relative importances of the various exchange diagrams depend on the supersymmet- 
ric parameters through the couplings and masses of the exchanged particles. Typically, one 
expects a small value of cr ann at small tan/3, where o ann is dominated by sfermion and Z 
exchanges. In fact, in supergravity models at small tan/3, Higgs-exchange contributions are 
reduced not only by small couplings, but also by large values of the mass Ma- These features 
are displayed in Fig. 2, which is for tan/3 = 1.5 and 8\ = 62 = 0. As expected, fi x h 2 is 
rather large, and many configurations are excluded by the cosmological bound Qh 2 < 1. The 
allowed region (denoted by squares) is mainly due to configurations where o ann is dominated 
by I exchange. Filled diamonds denote configurations where Q x h 2 is in the range of Eq.(|j). 
In these configurations typical values for the mass parameters are: mj ~ (150 — 200) GeV, 
m q ~ 600 GeV, Ma — 1 TeV. In Fig. 2(b) we notice that the allowed domain extends to 
the right of the line m = m 0tmin (where m 0tmin is the minimal value for m ), allowing for 
the neutralino only a gaugino-dominated region. This occurs because the coefficient J2 of 
Eq.(^) is positive (see Fig.l). 

An illustration of how a deviation from mo universality may somewhat modify the picture 
is provided by Fig. 3, where we display the relevant physical regions when we choose a 
departure from universality: 5i = —1.0,52 = 1-0, which makes the coefficient J2 very small 
and negative. As a consequence, m^j and [i are strongly correlated, with a slight extension of 
the neutralino physical region toward the sector of higgsino-gaugino mixture for m 1( / 2 > 400 
GeV (see Fig. 3(b)). These new configurations, which turn out to correspond to high values 
of mo: mo — (1.3 — 1.6) TeV, are no longer in conflict with the cosmological bound, at 
variance with the universal case. In these points of the parameter space a ann is dominated 
by the i-exchange, since here m- t ~ 250 GeV. However, we remark that these configurations 
realize radiative EWSB only with strong fine tuning. The generic trend is to have, for 
5\ = —1.0, 62 = 1.0, Q x h 2 smaller than in the universal case and this feature is displayed in 
Fig. 3, where we notice that the cosmological constraint is less effective here in restricting 
the parameter space, as compared to the universal case (see Fig. 2). 



Let us turn now to large values of tan/3. As was already pointed out in Ref. [14] for 



the universal case at tan/? = 53, the neutralino relic abundance is very low: Q x h 2 < 0.1, 
due to small values of Ma- Ma < 150 GeV, which generates a large a ann dominated by 
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Higgs-exchange contributions. 

By moving away from tuq universality, we can change the picture noticeably. For instance, 
by taking 8\ = 0.4,^2 = —0.1, we can generate large values of Q x h 2 . However, as we see 
in Fig. 4, most of the configurations of large Q x h 2 are disallowed by the constraint on the 
bottom mass (and also marginally by the constraint b — > 57). The reason for the marked 
difference between this case and the universal one is that now M4 is large. Indeed, the 
lowest allowed values of Ma {Ma — 150 GeV) only occur in the parameter region: mo — 
300 GeV, mi/2 — 400 GeV. The large values of M4 are caused by a sizeable, positive value 
for the coefficient K 2 of Eq.(|7]). The only allowed configurations which provide Q x h 2 in the 
range of Eq.(|) are those with m 1/2 ~ 120 GeV and 1.2 TeV < m < 1.8 TeV. Notice that 
some of these configurations do not satisfy our no- fine-tuning constraint. The relevance of 
the a s value in determining Q x h 2 and in shaping the allowed physical regions is apparent 
when Fig. 4, evaluated for a s = 0.1127, is compared to Fig. 5, evaluated for ot s = 0.118. 

V. EVALUATION OF THE SIGNALS 

As was anticipated in the Introduction, the signals to be discussed in the present paper 
consist of the fluxes of up-going muons through a neutrino telescope generated by neutrinos 
which are produced by pair annihilations of neutralinos captured and accumulated inside 
the Earth and the Sun. The steps involved in this process are the following: a) capture 
by the celestial body of the relic neutralinos through a slow-down process due essentially 
to neutralino elastic scattering off the nuclei of the macroscopic body, b) accumulation of 
the captured neutralinos in the central part of the celestial body, c) neutralino-neutralino 
annihilation with emission of neutrinos, and for the various annihilation products, d) propa- 
gation of neutrinos and conversion of their component into muons in the rock surrounding 
the detector (or, much less efficiently, inside the detector), and finally e) propagation and 
detection of the ensuing up-going muons in the detector. 

The various quantities relevant for the previous steps are calculated here according to 
the method described in Ref. ||, to which we refer for further details. 

A. Neutralino local density 

As will be shown explicitly in the subsequent subsection, the capture rate of the relic 
neutralinos is proportional to the local neutralino density p x in the solar neighbourhood. 
Let us specify here how p x is evaluated. 

For each point of the model parameter space we first calculate the relevant value of 
the cosmo logical neutralino relic density according to the procedure outlined in Sect. 4. 
Whenever Q x h 2 is larger than a minimal (Qh 2 ) min suggested by observational data and by 
large scale structure calculations, we simply put p x = pi, where pi is the local halo density. In 
the points of the parameter space where Q x h 2 is less than (Qh 2 ) min , the neutralino may only 
provide a fractional contribution Q x h 2 / (Qh 2 ) m i n = £ to Qh 2 ; in this case we take p x = pi^. 
The value to be assigned to (Qh 2 ) m i n is somewhat arbitrary. Here it is set equal to 0.1. As far 
as the local halo density pi is concerned, we are inspired by the recent estimate pi = 0.51 jZo!i7 
GeV cm -3 , based on a flattened dark matter distribution and recent microlensing data. 
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This introduces a significantly larger central value as compared to previous determinations 
(see, for instance, Ref. [0]). All the numerical results presented in this paper were obtained 
using the value p\ = 0.5 GeV cm -3 . 



B. Capture rates and annihilation rates 

For the evaluation of the capture rate C of the relic neutralinos by a celestial body we 



have used the standard formula 38 



c = rE^rWi)(«D^ (8) 



v x i m x m 



where v x is the neutralino mean velocity, <7j is the cross section of the neutralino elastic 
scattering off the nucleus % of mass m„ MbU is the total mass of the element % in the 
body of mass Mb, (vl sc )i is the square escape velocity averaged over the distribution of the 
element i, and X, is a factor which takes account of kinematical properties occurring in the 
neutralino-nucleus interactions. For the evaluation of the elastic x~ nuc l eus cross sections 



we refer to 14 



The annihilation rate Ta of the neutralinos inside the macroscopic body is calculated 



with the formula 39 



C , 2 ft 



T A = -tanh 2 — (9) 

I \T A ) 

where t is the age of the macroscopic body [t = 4.5 Gyr for Sun and Earth), ta = (CCa)^ 1 ^ 2 , 
and Ca is the annihilation rate per effective volume of the body, given by 

ca - ^ (M 3 ' 2 (10) 



V V20 GeV 

Here, V is defined as V = (3m 2 Pl T/(2p x 10 GeV)) 3/2 where T and p are the central 
temperature and the central density of the celestial body. For the Earth (T = 6000 K, 
p = 13 g • cm" 3 ) V = 2.3 x 10 25 cm 3 , for the Sun (T = 1.4 x 10 7 K, p = 150 g • cm- 3 ) 
Vo = 6.6 x 10 28 cm 3 . Also, <7 ann is the neutralino-neutralino annihilation cross section and v 
is the relative velocity: < er ann t> > is calculated with all the contributions at the tree level as 
previously discussed in Sect. 4, with the further inclusion here of the two-gluon annihilation 
final state [HO . 



From Eq.(||) it follows that in a given macroscopic body the equilibrium between capture 
and annihilation (i.e. Ta ~ C/2 ) is established only when t > ta- We stress here that the 
neutralino density p x enters not only in C but also in ta (through C). Thus, the use of a 
correct value for p x (rescaled, when necessary) is important also in determining whether or 
not equilibrium is reached in a macroscopic body. 

From the evaluation of the annihilation rate for neutralinos inside the Earth and the 
Sun it turns out that, for the Earth, the equilibrium condition depends sensitively on the 
values of the model parameters, and is not satisfied in wide regions of the parameter space. 
Consequently, for these regions the signal due to neutralino annihilation may be significantly 
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attenuated. On the contrary, in the case of the Sun, equilibrium between capture and 
annihilation is reached for the whole range of m x , due to the much more efficient capture 
rate due to the stronger gravitational field p8 |. 



C. Neutrino fluxes 

Let us turn now to the evaluation of the neutrino fluxes due to the annihilation processes 
taking place in the celestial bodies. For a distant source such as the Sun the differential rate 
in the neutrino energy E u is given by 

dN v = T A B (F) dN fu 

dE v 4ird 2 frj xf dE v 1 ; 

( F) 

where d is the distance from the source, F denotes the x~X annihilation final states, B^J 
denotes the branching ratios into heavy quarks, r leptons and gluons in the channel F; 
dNf u /dE u is the differential distribution of the neutrinos generated by the hadronization of 
quarks and gluons and the subsequent hadronic semileptonic decays. 

In the case of the Earth one has to take into account the size of the region around the 
center of the Earth where most of the neutralinos are accumulated. This is important, since 
the angular dependence of the flux plays a crucial role in providing a signature, and, po- 
tentially, in allowing a determination of the neutralino mass [08j|6i^T|. For details about the 
relevant formulae we refer to ||. Here we only give some information about our evaluation 
of the spectrum dNf u /dE u to be employed in Eq. flTTD . The neutrino spectra due to b and c 
quarks, r leptons and gluons were computed using the Jetset 7.2 Monte Carlo code We 



have neglected the contributions of the light quarks directly produced in the annihilation 
process or in the hadronization of heavy quarks and gluons, because these light particles 
stop inside the medium (Sun or Earth) before their decay ||43|| . For the case of the Sun we 



have also considered the energy loss of the heavy hadrons in the solar medium. The spectra 
due to heavier final states, i.e. Higgs bosons, gauge bosons and t quark, were computed 
analytically by following the decay chain down to the production of a b quark, c quark or 
a tau lepton; the result of the Monte Carlo was used to obtain the final neutrino output. 
Because of the high column density of the solar medium, the absorption and the energy loss 
of the produced neutrinos were also included. 

One possible effect that we have not included is that of matter-enhanced neutrino os- 
cillations. As pointed out in Ref. MM, these could be important for neutrinos produced by 



neutralino annihilations in the Sun, if there is a large mixing angle as in some fits to solar 
or atmospheric neutrino data. It would have been more conservative to allow for z/ M — > v e 
or v T oscillations inside the Sun, but we neglect them here. 



D. Fluxes of up-going muons 

The capability of a neutrino telescope to measure the flux of Eq.flTTD depends on how 
well this signal may be discriminated from the background due to neutrinos produced in the 
Earth's atmosphere by cosmic rays. In the case of the Sun, the signal to background (S/B) 
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discrimination is based on the correlation with the position of the Sun in the sky. As far as 
the signal from the Earth is concerned, the S/B discrimination requires an analysis of the 
angular distribution, and is based on the property that, whereas the angular distribution 
of the signal is expected to be markedly peaked at the nadir, the background distribution 
is almost flat (with a slight increase at the horizon). In order to evaluate these angular 
distributions properly, one has to take into account the size of the region around the center 
of the Earth where neutralinos are expected to accumulate, and the physical processes that 
entail an angular spreading of the signal. Some angular spreading is induced by the — > \i 
conversion: 9^ v ~ 2°(^/100 GeV)- 1 / 2 , and a comparable effect occurs because of the 
multiple scattering of the muon in the rock surrounding the detector. Thus, an accurate 
calculation of angular distributions requires a Monte Carlo simulation ||. 

In this paper, since we only deal with spectra, and not with angular distributions, we 
perform the calculation of the up-going muon fluxes in a simple no-straggling approximation 
In fact, this approximation reproduces accurately the Monte Carlo calculations for the 



distributions in the muon energy H. Then, we can write the muon spectrum as [45 



dE; ~ Na aVb% k >th dK iK k " dEi (12) 

where E^ th is the muon energy threshold, is the Avogadro number, da(E u , E' ) /dE' is 
the differential cross section for the production of a /i from a impinging on an isoscalar 
target and A + BE^ is the average muon energy loss, due to ionization, pair production, 
bremsstrahlung and photonuclear effects. For the coefficients A and B we have used the 
following values A = 2.4 • 1(T 3 GeV/(g/cm 2 ), B = 4.75 • KT 6 (g/cm 2 )" 1 [||. 



Let us conclude this section with some comments on the typical parameters of the re- 
quired experimental layout. In deep underwater /ice experiments muons from neutralino 
annihilation are planned to be detected using their direct Cerenkov radiation. A relativistic 
muon in water emits per 1 cm about 250 Cerenkov photons with wavelengths in the interval 
300 — 500 nm. Then it is easy to estimate the number of photoelectrons produced in a 
photomultipler (PM) at a distance r from muon trajectory as 

N e ~ 30^D 2 M (13) 
r 

where Q ~ 0.15—0.25 is the PM quantum efficiency, all distances are given in cm and we take 
Dpm — 35 cm as the PM diameter. One can see from Eq. (p~3| ) that at a distance equal to the 
scattering length in the water/ice, which we take as l sc ~ 20 — 30 m, PM's detect a strong 
signal corresponding to about 3 photoelectrons. To determine the muon trajectory, only the 
PM's located closer than the scattering length can be used. Their number, Npm — l^lc/d 3 , 
is of order 100(^/(100 GeV)) for a distance between detectors d — 10 m and a scattering 
length 25 m. 

The supersymmetric model we are discussing is characterized by light neutralinos, m x < 
200 GeV. This results in a small pathlength of the produced muons, 1^ = 500 (J5 M / (100 GeV)) 
m, which has the following observational consequences. 

• For a widely-discussed 1 Km 3 array, many muons have trajectories confined within the 
detector. The energies of these muons can be estimated from the trajectory lengths. 
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• There is a small probability to find a big shower due to bremsstrahlung or pair pro- 
duction along the muon trajectory. 

• A reliable measurement of the muon trajectory, which is necessary for identification 
of a neutralino-produced muon, needs a rather dense array with a distance between 
detectors d ~ 10 m in the case of a scattering length 25 - 30 m. It could be the core 
of a larger detector. 



VI. RESULTS AND CONCLUSIONS 

In this section we show our results for integrated fluxes $ M of the up-going muons, and 
compare them with the upper bounds obtained at the present neutrino telescopes. We will 
then compare these indirect signals with the signals measurable by direct searches for relic 
neutralinos. We also show how the strengths of the various signals correlate with each other 
and with the values of the neutralino relic abundance. However, we wish first to comment 
on some general properties which later provide a simple interpretation of some characteristic 
features of our numerical results. 

(i) Both the event rate for neutralino direct detection R and the capture rate C of a relic 
neutralino by a celestial body are proportional to p x cxj. This property also applies to $ /t , with 
some further more complicated dependence on Oi and cr ann for the flux $ M from the Earth, 
when equilibrium is not yet established (see Sect. 5.B). Thus, we expect R and $ M to show 
some similarities in their behaviours as functions of the supersymmetric parameters, when 
the same elastic cross sections <Tj's are involved. This is roughly the case, when we consider 
the indirect signal from the Earth. In fact, in this instance the neutralino capture from the 
macroscopic body occurs dominantly (except for extremely pure gaugino compositions of the 
neutralino) through coherent cross sections, which have an overwhelming role also in direct 
detection by the nuclei considered here. Obviously, this is a qualitative argument, since the 
difference in the nuclear compositions of the detectors and of the Earth play an important 
role in the actual numerical results. Furthermore, the similarity between R and ^ arth is 
attenuated when equilibrium in the Earth is not yet established, as previously mentioned. 
A correlation between R and $^ arih is manifest in some of the figures that will be presented 
below. At variance with the case of the Earth, it has to be noted that, in some regions of 
the supersymmetric parameter space, neutralino capture by the Sun occurs mainly through 
spin-dependent cross sections, due to the overwhelming presence of Hydrogen in the Sun. 
Under these circumstances, R and the signal from the Sun $^ nn are not expected necessarily 
to resemble each other. However, signals which are not enhanced by coherent effects, are 
generally much below the experimental sensitivities. Summarizing, it is worth remarking 
that, in view of the previous arguments, direct detection and indirect detection by up-going 
muons have some common regions of the parameter space to explore, but may also have 
very distinctive features in their discovery potentials |I7 



(ii) If we combine the properties: R oc pi<Ji, ^ oc pidi, and Vt x h 2 oc< a^v >^J, and take 
into account the fact that usually <7j and cr ann , as functions of the supersymmetric model 
parameters, are either both increasing or both decreasing, we come to the conclusion that the 
quantities $ M and R are somewhat anticorrelated with £l x h 2 (see also P%|]). This property 
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is attenuated, but usually not washed out, in the case of scaling of the neutralino local 
density. In fact, when scaling occurs (i.e., when fl x h 2 < (Qh 2 ) min ), one has R oc p/£o"j oc 
Pi<Ji x [< cr ann t> >int]~\ an d the same for Since it turns out that when cr, is large also <r ann 
increases, but in such a way that usually the ratio crj/a ann increases too, one can conclude 
that a form of anticorrelation between fl x h 2 and ^ (or R) persists also when scaling of p\ 
is effective. A feature of this type is displayed by our numerical results. 

Let us turn now to a presentation of some of our results. We show in Fig. 6 a sample 
of our calculations for the signals expected from the Earth at the representative point tan/? 
= 53, for three different choices of the 5i parameters. The three closed curves denote the 
boundaries of the allowed regions when the parameters mo and mi/2 are varied in the ranges 
10 GeV < mo < 2 TeV, 45 GeV < mu2 < 500 GeV. The figure displays the flux of up-going 
muons integrated over a cone of half aperture of 30° centered at the nadir, and for muon 
energies above 1 GeV. Also shown in Fig. 6 is the most stringent experimental bound: 
^Earth < 2.1 • 10" 14 cm" 2 - s- 1 (90% C.L.), obtained with an exposure of 2954 m 2 - yr fT§ 
(for other recent experimental data see [^9|-|5T[|). 

In Fig. 6(a) we notice, in the case d~i = 0, 62 = 0, a s = 0.1127 indicated by the solid 
line, a decreasing behaviour of $ M as a function of the neutralino mass, as expected from 
the structure of the capture rate C. It is remarkable that for some neutralino configurations 
in the interesting range m x ~ 150 GeV the present experimental sensitivity (or a slight 
improvement of it) is already adequate for restricting the supersymmetric parameter space. 
These configurations have a light A: Ma < 70 GeV. The correlations between $ M and £l x h 2 
or between $ M and R, anticipated in points (i) and (ii) above, are apparent in parts (b) 
and (c) of Fig. 6, respectively. One important feature of this parameter choice is shown 
by the solid line in Fig. 6(b): the neutralino configurations, which are explorable with 
measurements of $ A1 , yield values of £l x h 2 below the range of Eq.(|J). We stress that this is a 
general trend, which will be further discussed at the end of this Section. Fig. 6(c) shows how 
measurements of $ M and R may give information about similar neutralino configurations. 
The event rate for direct detection R, employed in the present paper as well as in our previous 



work [ 141 , is defined as the integral of dR/dE ee (E ee is the electron equivalent energy) over 
the 12-13 KeV range, for a Ge detector. The experimental upper bound is that obtained 
from the experiment of Ref . [|52| (for other experiments using Ge detectors see Ref. ]53| ) . 

The discovery potential of neutrino telescopes in this context is further illustrated in 
Fig. 7, where the case denoted by the solid line in Fig 6 is shown. This figure displays (by 
empty squares) the region in the m 1( / 2 — m plot (and in the mi/ 2 — p plot) which may be 
explored by a neutrino telescope whose sensitivity is one order of magnitude better than that 
of Ref. |n| . The region denoted by light crosses would require an even better sensitivity. In 
this connection, we stress that the theoretical results presented here were obtained using a 
particular set of representative values, not only for some supersymmetric model parameters, 
but also for a number of astrophysical and cosmological parameters, many of which are 
affected by large experimental uncertainties. This is the case, for instance, of the dark 
matter local density pi and of Qh 2 . For this reason, one has to be extremely cautious in 
extracting exclusion plots in the supersymmetric parameters from a comparison between 
present experimental and theoretical values of $ At . Obviously, a region of the neutralino 
parameter space would be excluded only if the theoretical values, obtained with the most 
conservative set of values for the free parameters, were above the experimental upper bound. 
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This is certainly not the case with the present experimental limits. Particular attention has 
to be paid to the value of the neutralino local density p x to be used. Just employing p x = pi 
everywhere, without making a consistency check with the value of Q x h 2 and without rescaling 
appropriately, leads to erroneous exclusion plots. 

In Fig. 6 we present as dashed lines another set of results, corresponding to a case where a 
departure from universality in soft scalar masses is introduced: 5\ = 0, 5 2 = —0.3. We notice 
that many configurations may produce signals above the experimental bounds, both in direct 
and indirect searches. The dashed lines exhibit correlations between the various quantities 
similar to those exhibited by the solid lines, although in this case some configurations provide 
both a measurable $ M and an Q x h 2 in the range of Eq. (|4]). The configurations whose 
signals are above the experimental limit are displayed in Fig. 8 with heavy oblique crosses. 
In particular, from Fig. 8(b) it is clear that the corresponding neutralino compositions are 
dominantly gaugino-like (notice that the allowed region in the m\j 2 -p plane extends to the 
right of the line m = m miri , due to a negative value of the coefficient J 2 of Eq.(Q)). It 
is worth recalling that the strict nii/ 2 -p correlation appearing in Fig. 7 is removed here, 
due to the deviation from universality we have introduced (the coefficient J 2 is positive and 
sizeable here). The reason why many configurations are above the experimental bound, in 
spite of their their gaugino nature, is the fact that in these configurations the CP-odd Higgs 
neutral boson A is light: < 65 GeV (the coefficient K 2 of Eq.(0) is negative). 

One further example of a non-universal case is provided by the dotted lines in Fig. 6 and 
by Fig. 9, where £1 = 0.7, 5 2 = 0.4. A number of configurations turn out again to exceed 
the experimental bound. Here the neutralino parameter space opens up to the left of the 
line m = mo^ min , where gaugino-higgsino mixture takes place. Again, the configurations 
with the highest signals have a light A boson. 

Examples of signals from the Sun, $^ nn , are shown in Figs. 10 and 11 together with the 
value of the experimental bound: $^ nn < 3.5 • 10~ 14 cm -2 - s _1 (90% C.L.), obtained with 
an exposure of 1002 m 2 - yr JTSj] . By comparing with Fig. 6, we notice that the distribution 
of points along the vertical axis is much more spread out in the case of the signal from the 
Earth than from the Sun, with many configurations giving values of ^^ arth far below (by 
many orders of magnitude) the experimental sensitivity. The strong attenuation of ^ arth 
for these configurations is due to the fact that for them the capture-annihilation equilibrium 
in the Earth is far from being established. As for the signals from the Sun, we further note 
that for all configurations of Fig. 10 the signal is dominated by coherent contributions, in 
spite of the scarcity of heavy nuclei in the Sun. This makes the discovery potential for the 
signal from the Sun essentially equivalent to the one for the signal from the Earth, at this 
representative point. Spin-dependent effects contribute only for a few per cent, at most. In 
other representative points, the main contribution to the signal may be due to non-coherent 
cross sections. However, under these circumstances the overall signal is much smaller than 
the experimental limit. 

In all the examples just discussed (Figs. 6-11), referring to the representative value 
tan/? = 53, we have found that measurements of $^ provides a useful tool for investigating 
interesting regions of the supersymmetric parameter space, when the boson A is light, Ma = 
O(Mz). Unfortunately, but unavoidably, the very fact that Ma is small makes it difficult for 
these supersymmetric configurations to yield a neutralino relic density in the desired range: 
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0.1 < Q x h 2 < 0.3. A quite opposite scenario occurs in different, and equally allowed, sectors 
of the parameter space, such as those examined in Sect. 4, where Q x h 2 is very large (and 
where the signals R and $ M are much below the present and the foreseeable experimental 
sensitivities). This is typically a scenario of small tan/5 or it may occur at large tan/5 with 
deviations from m universality (see, for instance, the case illustrated in Fig. 3). It is 
remarkable that both scenarios are possible in the case of the supergravity-inspired model 
adopted in the present paper, although it is constrained by the requirement of radiative 
EWSB. Of course, much wider regions of the parameter space are allowed to both scenarios 
in the case of the unconstrained MSSM. 
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Figure Captions 



Figure 1 - Coefficients J\, J2 and K\, K% of the polynomial expressions (6,7) as functions 
of tan /3. The range 1.5 < tan (3 < 4 is shown in section (a) and the range 4 < tan (3 < 53 in 
section (b). In the upper part Ji is denoted by a solid line, J2 by a dot-dashed line. In the 
lower part, K\ is denoted by a solid line, K 2 by a dot-dashed line. 



Figure 2 - a) The parameter space in the (m^, mo) plane for tan/? = 1.5, 5\ = 
and 62 = 0, a s = 0.1127. The empty regions are excluded by: i) accelerator constraints, ii) 
radiative EWSB conditions, hi) the LSP not being a neutralino. Dots represent the region 
where Q x h 2 > 1. Regions with crosses are excluded by b —>■ sj and constraints. In the 
regions denoted by filled diamonds, 0.1 < Q x h 2 < 0.3. b) The parameter space represented 
in the (m^, fi) plane. The line mo = mo )mm corresponds to the minimum value of mo. 
Notations are the same as in a), but crosses are omitted here. The dashed lines denote the 
no-fine-tuning upper bounds on m , mi/2 and fi. 

Figure 3 - The same as in Figure 2, but with 5\ = —1.0 and 62 = 1.0. 



Figure 4 - The same as in Figure 2, but with tan/3 = 53, 5\ = 0.4 and 62 = —0.1. 
a, = 0.1127. 



Figure 5 - The same as in Figure 4, but with a = 0.118. 



Figure 6 - The flux ^ arth for tan/3 = 53: 8 X = 0, 8 2 = and a s = 0.1127 (solid line), 
Si = 0, 62 = —0.3 and a s = 0.118 (dashed line), and 61 = 0.7, 62 = 0.4 and a s = 0.115 
(dotted line) as a function of m x (a), as a function of Q x h 2 (b), and plotted versus the rate 
for direct detection with a Ge detector (c). Allowed configurations stay inside the closed 
curves. The horizontal line displays the experimental bound ^ aTth < 2.1 ■ 10 -14 cm" 2 - s _1 
(90% C.L.) fll5| , and the vertical line of section (c) displays the upper limit of Ref. fl52 



Parameters are varied on a linear equally-spaced grid over the ranges: 10 GeV < mo < 
2 TeV, 45 GeV < m 1/2 < 500 GeV. 



Figure 7 - a) The parameter space in the (m^, m ) plane for tan/3 = 53, 5i = and 
62 = 0, a s = 0.1127. Empty regions are excluded by: i) accelerator constraints, ii) radiative 
EWSB conditions, iii) the LSP not being a neutralino, iv) b — > 37 and rrib constraints. 



The heavy oblique crosses denote configurations above the present experimental bound [15 



Squares denote regions which could be explored by a neutrino telescope whose sensitivity is 
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one order of magnitude better than that of Ref. [15]. The region with light crosses would 
require an even better sensitivity, b) The parameter space represented in the (m^, fi) 
plane. 



Figure 8 - The same as in Figure 7, but with 8\ = and 82 = —0.3, a s = 0.118. 



Figure 9 - The same as in Figure 8, but with 8\ = 0.7 and 82 = 0.4, a s = 0.115. 



Figure 10 - The flux $f un for tan f3 = 53, 8x = and 8 2 = -0.3, a = 0.118, as a function 



"I" 



of m x (a), as a function of Q x h 2 (b), and plotted versus the rate for direct detection with 
a Ge detector (c). Allowed configurations stay inside the closed curves. The horizontal line 
displays the experimental bound $^ un < 3.5-10 -14 cm" 2 - s _1 (90% C.L.) [|15|| , and the vertical 
line of section (c) displays the upper limit of Ref. [|52 . Parameters are varied on a linear 



equally-spaced grid over the ranges: 10 GeV < m < 2 TeV, 45 GeV < m^j < 500 GeV. 



Figure 11 - a) The parameter space in the (777.1/2, mo) plane for tan/3 = 53, 81 = 
and 82 = —0.3, a s = 0.118. Empty regions are excluded by: i) accelerator constraints, 
ii) radiative EWSB conditions, iii) the LSP not being a neutralino, iv) b —>■ 37 and nib 
constraints. The heavy oblique crosses denote configurations above the present experimental 
bound: $^ un < 3.5 • 10~ 14 cm -2 - s~ x (90% C.L.) | 15| . Squares denote regions which could 
be explored by a neutrino telescope whose sensitivity is one order of magnitude better than 
that of Ref. |TS[ . The region with light crosses would require an even better sensitivity, b) 



The parameter space represented in the (mi/2, /i) plane. 
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